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ABSTRACT 



The general unloading phase of amphibious operations is 
examined with particular emphasis on the queueing problems that 
arise. A model is structured and a computer simulation is pro- 
vided. Variables considered are the number of transfer vehicles 
and their speed and payload, the number and capacity of loading 
and unloading points, and the ship- to- shore distance. Basic 
relationships between the variables are examined and an approxi- 
mating formula is developed for computing the , time required to 
unload a given tonnage of cargo. Applications in operational 
planning, system modification, and total system design are out- 
lined. The impact on system performance of various assumed 

\ 

distributions of travel time and loading time is discussed. 

For systems employing mixed fleets of transfer vehicles, some 
basic decision rules are developed as to which vehicles to load 
or unload first. Suggestions as to further research and as to 
the model *s general applicability are included. 
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CHAPTER I 



INTRODUCTION AND RATIONALE 

One of the key characteristics that serve to distinguish 
amphibious operations from other types of military operations is the 
requirement to project forces ashore onto a hostile land mass. This 
necessitates the rapid and efficient landing of large numbers of men 
and large quantities of supplies and equipment. Two quite distinct 
problems are involved. The first is the tactical landing itself, 
during which the scheduled waves and the on-call serials are landed. 

The scheduled waves consist of troops and equipment boated in landing 
craft, amphibians, or helicopters. They are organized into waves and 
land according to a predetermined time schedule. The on-call serials 
consist of groupings of troops, equipment, or supplies, for which an 
early requirement ashore is anticipated. They are called ashore as 
the need for them arises. 

After the scheduled waves and the serials have been landed, 
the general unloading commences. During the general unloading phase, 
all remaining cargo and equipment is unloaded from the ships. Landing 
craft, amphibians, and helicopters, which we shall refer to as trans- 
fer vehicles, shuttle back and forth between the ships and the beach- 
head, bringing ashore the large quantities of supplies necessary to 
sustain the combat forces which have been landed. This is not a glam- 
orous or exciting phase of an amphibious operation, but the importance 
of a rapid and efficient build-up of supplies ashore to maintain the 
momentum of the attack cannot be overlooked. 

As the transfer vehicles shuttle back and forth, they may some- 
times have to wait at the ship while other vehicles are being loaded. 
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Likewise, they may have to wait at the beach while other transfer 
vehicles are being unloaded. The entire system is circular in nature, 
and may be schematically represented as in Figure 1. In making a cycle, 
each transfer vehicle may go through six successive stages; the loading 
queue, the loading operation itself, a transit to the unloading area, 
the unloading queue, the unloading operation, and a transit back to the 
loading area. Of course, sometimes an arriving vehicle will find that 
there is no queue of vehicles waiting to load or unload, as the case may 
be, and it will move directly into the loading or the unloading operation. 

Many factors influence the way in which this system operates, 
including transfer vehicle speed and payload, the ship-to-shore distance, 
loading and unloading rates, and the number of transfer vehicles in the 
system. In analyzing the effects of these various factors, it becomes 
apparent that we cannot neglect the queueing problems which are inherent 
in the physical system being considered. Unfortunately, the system which 
this paper will model does not lend itself to solution by the usual queue- 
^ ing theory methods. The usual queueing theory models assume that arrivals 
into the system come from some infinite population at a rate which can 
be measured and in accordance with a probability distribution which can 
be determined. The system we wish to model, however, involves a finite 
number of transfer vehicles. Their number may, in fact, be quite small 
in many actual operations. The rate at which the transfer vehicles 
arrive at the loading or unloading areas for service is dependent on 
what is happening at all other stages in the system. Additional compli- 
cations will arise if we wish to consider the transient period which 
exists before the system settles down to a steady state. A particularly 
difficult problem arises from the fact that in practice the fleet of 
transfer vehicles will often not consist of a single type of vehicle. 
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FIGURE 1 

Schematic Representation of 
General Unloading Operation 
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but rather will consist of a variety of vehicles, perhaps with widely 
different operating characteristics. 

In this paper we will structure the queueing system and describe 

and discuss the results of a simulation through which the relationships 

between several of the variables in this system have been analyzed, and 

by means of which some of the theoretical questions concerning the nature 

of this queueing system have been investigated. In a recent paper, 

12 

Posner presented a closed-form method for computing the steady-state 
probabilities for queueing systems of the type which we wish to model, 
However, the computations necessary in his approach are lengthy and 
solutions are dependent upon the assumption of exponential service rates. 
A computer simulation, on the other hand, will remove the need for this 
assumption and also enables us to examine the effect of the transient 
period at the start of system operation. The computer simulation also 
permits consideration of various priorities of service when dealing with 
a fleet of transfer vehicles consisting of a variety of vehicle types. 

Chapter II will present a description of the model and the com- 
puter simulation will be described in Chapter III, In Chapter IV, 
results from a series of simulation runs will be presented and analyzed. 
The parameter values which are used do not correspond to any actual 
equipment, but are representative of a range of values attainable by 
existing or proposed equipment. The results serve to demonstrate the 
type of information that may be generated by the simulation and the 
general nature of the relationships between the various parameters. 

An empirically derived formula for estimating the total time to com- 
plete the unloading of a given tonnage of cargo will be presented. In 
Chapter V, the simulation program will be used to test whether the sys- 
tem exhibits the characteristics of a Poisson process, under various 
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assumptions about the underlying distributions of travel time and 
service time. An hypothesized mean arrival rate at the loading area 
will also be tested. This will give us some insight into the nature 
of the system being studied and will provide us with information as 
to the significance of the underlying probability distributions on 
overall system performance. Some results obtained employing mixed 
fleets of transfer vehicles will be presented in Chapter VI. A sum- 
mary of results thus far obtained, conclusions, and some comments as 
to the applicability of the model to other systems completes this paper. 
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CHAPTER II 



DESCRIPTION OF THE MODEL 

The physical system we wish to model is the general unloading 
phase of an amphibious operation. The system may be characterized as 
containing several key elements. First, there are the ocean-going 
ships which have transported the landing forces and their supplies 
and equipment to the objective area. We assume that port facilites 
are not available for direct unloading onto piers. Next, we have a 
fleet of transfer vehicles to move the cargo from the ships to the 
shore. These transfer vehicles may be landing craft, amphibians, heli- 
copters, surface-effect machines, or hydrofoils. In actual operations 
they will probably exist as a mixed fleet of various types. Whatever 
their physical characteristics, however, they can be categorized accord- 
ing to their speed and their payload. In this model, we define speed 
to be a readily attainable cruising speed, rather than the vehicle’s 
top speed. Thus, we may meaningfully speak of a mean time to travel 
a given distance. We assume that the vehicle’s speed is the same, 
whether it is loaded or unloaded. Payload is taken to be a usable pay- 
load, readily achievable when the vehicle is loaded with standardized 
cargo, such as loaded pallets of uniform size and weight. 

Other elements in the system are the cargo-handling facilities 
on the ships and on the shore. These facilities may also have widely 
differing physical characteristics, but they can be categorized by 
their capacity and by the number of loading points. By capacity, we 
mean their cargo-transfer rate, in tons per hour, achieved when load- 
ing standardized cargo into the type of vehicle under consideration. 

A loading point is defined as a location and a cargo-handling facility 
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that can accommodate one transfer vehicle at a time. Two possible 
queues in the system are a queue of transfer vehicles waiting at the 
ship to load cargo and a queue of transfer vehicles at the unloading 
area waiting to unload. 

Characteristics of the Model . 

A complete model of the general unloading operation should 
include all the steps involved in moving cargo from the hold of a ship 
to an inland dump, or perhaps all the way to the user. However, the 
model developed for this study has been simplified in that it considers 
only a movement of cargo from the ship to an unloading area, which may 
be thought of as the waterline at the beach, although in the case of 
amphibian vehicles or helicopters, it might be an inland dump. In 
either case, we shall assume that the transfer vehicles do not encounter 
any traffic delays enroute between the ship and the unloading area. The 
model also assumes that the vehicles cycle continually through the sys- 
tem until all the cargo has been moved. In other words, no allowance 
is made for refueling or maintenance time during the operation. 

In the computation of the loading time or the unloading time 
for a vehicle, no allowance is presently made in this study for the 
buffer time, or the time to move a vehicle into position. Buffer time, 
however, could be taken into account if desired by adding a constant to 
the computed loading or unloading time. Our computation of loading 
(unloading) time assumes that the time is a linear function of vehicle 
payload. This seems reasonable if the cargo to be moved consists of 
standardized cargo, since the time will be determined largely by the 
number of slingloads or cycles that the cargo-handling equipment will 
have to perform. However, this method of computing loading and unload- 
ing times may not be completely accurate, since the larger payload 
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vehicles may be easier to load or unload because of larger loading 
area, wider hatches, or other physical characteristics. In any com- 
parison of actual vehicles, if one vehicle lends itself to more expe- 
ditious loading, then this fact should be taken into account when 
computing loading or unloading time. 

A basic result from queueing theory is the fact that a single 
queue feeding customers (in this case transfer vehicles) into multiple 

i 

servicing facilities (in this case the loading or unloading points) 
will be more efficient than having a separate queue for each servicing 
facility. Therefore, the model is structured in this fashion, and no 
provision is made for other arrangements of the queues. 

Independent and Dependent Variables . 

Among the many variables which affect the system being con- 
sidered by this model, those which are of principal interest are 

1) the number of transfer vehicles, 

2) the’ speed of the transfer vehicles, 

3) the payload of the transfer vehicles, 

4) the distance the cargo must be moved, 

5) the number of loading and unloading points, and 

6) the rate at which the cargo-handling facilities can 
. load or unload the cargo. 

The model and the associated computer simulation are designed 
to treat all of these factors as independent variables, that is, values 
for these variables must be assigned for the program. The principal 
dependent -variables , whose values are the output of the simulation, are 

1) the time to complete the unloading, 

2) the avi'.rage^ waiting times per vehicle round trip. 
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3) the probabilities that the queues are of various 
lengths, and 

4) the amount of time the loading and unloading facili- 
ties are not completely utilized. 

As mentioned above, the cargo itself is considered to be 
standardized cargo consisting of pallets of uniform size and weight, 
and the model is not designed to consider the effects of variations 
in the cargo itself. This assumption about the uniformity of the 
cargo does not seem unreasonable, since the actual cargo involved in 
general unloading operations consists largely of rations, ammunition, 
and boxed military equipment, all of which is palletized in a standard 
and fairly uniform manner. 

Measures of Effectiveness . 

In this study, the time to complete the unloading of a fixed 
tonnage of cargo is taken as the basic measure of effectiveness. Any 
system that takes less time than another system is considered ''better** 
than that other system. Many different combinations of parameter 
values will produce approximately the same overall results, in terms 
of time to move a fixed tonnage of cargo. Changing the value of any 
one of the parameters will change the total time required to complete 
the operation. In any given situation, total time may be very sensi- 
tive to changes in one or more of the parameters and relatively insensi- 
tive to changes in other parameters. Trade-offs exist among the six 
independent variables which we are considering. 

In a study of broader range, of course, we would be interested 
not only in the time to complete the unloading operation, but also in 
the cost. However, it will not be possible to simultaneously minimize 
time and cost, since it appears that a higher cost system (say, a system 
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with more transfer vehicles) will always achieve at least some reduc- 
tion in time. Therefore, it will be necessary to minimize the cost 
to complete the unloading in a given time, or to minimize the time 
to complete the operation for a fixed cost. 

Using time as the basic measure of effectiveness seems intu- 
itively appealing, since time is certainly of the essence in an amphib- 
ious operation. Other possible measures of effectiveness, such as 
average queueing time or utilization rate of the cargo-handling facil- 
ities, do not appear to be of central importance. At any rate, their 
impact on system performance could be taken into account by extending 
the study to the costing of alternative systems. 
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CHAPTER III 



DESCRIPTION OF THE COMPUTER SIMULATION 

The computer simulation written for this study is and event- 
store simulation written in FORTRAN IV for the IBM 360/67. Program 
documentation is included in the appendices. Throughout the program, 
the time required to complete an event is computed in the following 
fashion. First, a random number is selected from a uniform distribu- 
tion, using the IBM library subroutine RANDU. The resulting uniform 
random number is then transformed to the desired probability distribu- 
tion. For example, in the case of a transformation from a uniform 
random number y to an exponentially distributed random number x, the 
formula is x=(-l/L)(ln y) , where 1/L is the mean of the exponential 
distribution. The number thus computed is then added to current time 
to obtain the time at which the event in question will occur. This 
time of occurrence is then stored in an array indexed by vehicle 
number . 

The program consists of seven major parts: 

1) the starting routine, 

2) the executive routine, 

3) arrival at the loading area, 

4) arrival at the unloading area, 

5) departure from the loading area, 

6) departure from the unloading area, and 

7) the ending routine. 

In the starting routine, constants and arrays are initialized and para- 
meter values are assigned to the independent variables. The simulation 
is started by placing all the transfer vehicles in the loading queue at 
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the ship, moving as many as possible into the loading points, and com- 
puting the time at which loading will be completed for each of the 
vehicles being loaded. 

The second part of the program, which corresponds to an execu- 
tive routine, searches through the arrays in which the occurrence times 
have been stored and selects the minimum time. This is the next event 
to occur, and can be one of four possibilities; an arrival at the load- 
ing area, an arrival at the unloading area, a departure from the load- 
ing area, or a departure from the unloading area. If the next event is 
an arrival at the loading area, the program checks the length of the 
queue and then the status of the loading points. If a queue exists, or 
if all loading points are occupied, the arriving vehicle goes into the 
queue. Otherwise it moves directly into an unoccupied loading point 
and a time for completion of loading is computed and stored. Arrivals 
at the unloading area are treated in the same fashion as arrivals at 
the loading area. 

If the next event is a departure from the loading area, a tran- 
sit time is computed and stored. The program then brings a vehicle 
from the queue into the loading point, or if no queue exists, sets the 
status of the loading point to empty. Departures from the unloading 
area are treated in the same fashion as departures from the loading 
area. Throughout the program, whenever queue length changes, the pro- 
gram computes the time the queue has been of a given length. 

When a specified quantity of cargo has been moved through the 
system, the program goes to its ending routine and computes and prints 
the data for that run. The output consists of the values of the depen- 
dent variables, as described earlier, together with identifying informa- 
tion for ease in analyzing the results. Modifications to the basic 
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program to accomplish specific tasks or to accumulate specific data 
will be discussed in the sections concerned with the special problems 
necessitating those changes. 
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CHAPTER IV 



DEVELOPMENT OF AN APPROXIMATING EQUATION 
FOR TIME REQUIRED TO UNLOAD A GIVEN TONNAGE OF CARGO 

As has been explained previously, the time required to complete 
the unloading of a fixed tonnage of cargo is the basic measure of effec- 
tiveness being used in this model. Six independent variables are being 
considered: the distance the cargo must be moved, the speed, the pay- 
load, and the number of transfer vehicles, the number of loading and 
unloading points, and the capacity of the loading and unloading points. 
Total tonnage of cargo to be moved is held constant at 10,000 tons 
throughout this paper, althougji the tonnage may easily be varied in 
the computer simulation. The simulation is structured as a balanced 
system, that is, with an equal number of loading and unloading points 
of equal capacity. Intuition, theory, and previous studies (ref. 1) 
all support the notion that an unbalanced system will not provide a 
desirable or "efficient** system. If, for instance, the total capacity 
of the loading points is greater than the total capacity of the unload- 
ing points, transfer vehicles will tend to queue up waiting to unload, 
thus producing inefficient utilization of the transfer vehicles. The 
loading points will be working at far less than their capacity, also, 
thus producing an inefficient overall system. 

In order to generate data for analysis, a number of computer 
runs were made with various values assigned to the independent vari- 
ables. The number of loading and unloading points was held constant 
at four each. Fleets consisting of eight, twelve, and sixteen transfer 
vehicles were considered. Four values of transfer vehicle speed were 
considered, in conjunction with four values of ship- to-shore distance. 
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as shown in Table I. In the computer simulation itself, only the mean 
travel time is used. Therefore we were able to consider all sixteen 
combinations of speed and distance by considering only seven mean 
travel times. 



TABLE I 

Mean Travel Times Utilized in Simulation Runs 

Ship-tO"Shore Distance 



Speed 




2 . 5 miles 


5 miles 


10 miles 


20 miles 


5 knots 




0 . 5 hrs • 


1.0 hr . 


2.0 hr. 


4.0 hr. 


10 knots 




0.25 hr. 


0.5 hr. 


1.0 hr. 


2.0 hr. 


20 knots 




0.125 hr. 


0.25 hr. 


0.5 hr. 


1.0 hr. 


40 knots 




0.0625 hr. 


0.125 hr. 


0.25 hr. 


0.5 hr. 



The mean capacity of each of the loading and unloading points was 
taken to be fifteen tons per hour, which appears to be in the range 
of the actual capacity of presently used equipment. Vehicle payloads 
considered were five, ten, fifteen, and twenty tons. For these runs, 
the underlying distribution of both service (loading or unloading) and 
travel times was assumed to be exponential. Chapter V discusses some 
results obtained using other assumed distributions. 

Portions of the sample data are presented graphically in 
Figures 2, 3, 4, and 5. In all of these figures, the principal depend- 
ent variable, time, is plotted on the vertical axis. This is the time 
to complete the movement of ten thousand tons of standardized cargo 
from the ships to the shore. In any graphical presentation of the 
data, only one or at most two of the independent variables can be dis- 
played conveniently in a single figure. In both Figures 2 and 3 the 
ship- to-shore distance is held constant at five miles, the number of 
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loading and unloading points is held constant at four each, and their 
mean capacity is held constant at fifteen tons per hour. In Figure 2, 
transfer vehicle speed is also held constant at ten knots and vehicle 
payload is plotted against time. Three separate curves are shown for 
various values of N, the number of transfer vehicles in the system. 

It is readily apparent from Figure 2 that, other things being equal, 
larger payload vehicles are more efficient than small payload vehicles. 
However, it is also apparent that the marginal returns from increasing 
vehicle payload diminish very rapidly. For a fleet of 12 vehicles, 
under the conditions portrayed by this figure, the marginal returns 
achieved by increasing vehicle payload from 5 tons to 10 tons are 
fairly substantial. Time to complete the operation is reduced from 
287 hours to 215 hours, for a gain of 72 hours. However, the marginal 
return achieved by increasing payload from 10 tons to 15 tons is only 
18 hours, and the marginal return achieved by increasing payload from 
15 tons to 20 tons is an insignificant 3 hours. 

In Figure 3, vehicle payload is held constant at 5 tons and 
vehicle speed is plotted against time. Again, three curves are shown, 
corresponding to 8, 12, and 16 vehicles in the system. Increasing 
vehicle speed from 5 knots to 10 knots produces impressive reductions 
in the total time necessary to move 10,000 tons of cargo through the 
system. Increasing vehicle speed from 10 to 20 knots produces a sub- 
stantial reduction in time for the fleet of 8 vehicles and lesser 
reductions for the fleets of 12 and 16 vehicles. As we increase 
vehicle speed above 20 knots, however, the marginal returns diminish 
almost to the vanishing point, as evidenced by the essentially hori- 
zontal shape of the curves between 20 knots and 40 knots. For the 
fleet of 16 vehicles, there is almost no time reduction achieved by 
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FIGURE 2 



Time to Unload Ten Thousand Tons of Cargo versus Transfer 
Vehicle Payload (Distance 5 Miles, Speed 10 Knots) 
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FIGURE 3 



Time to Unload Ten Thousand Tons of Cargo versus Transfer 
Vehicle Speed (Distance 5 Miles, Payload 5 Tons) 
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increasing speed from 20 to 40 knots under the set of conditions por- 
trayed in this figure. We would normally expect to encounter dimin- 
ishing marginal returns whenever we increase one parameter while 
holding the other parameters constant. In this system, however, the 
marginal returns from increasing vehicle payload, speed, or the num- 
ber of vehicles diminish very rapidly once the queue lengths begin 
to increase. Conversely, marginal returns from increasing the number 
or the capacity of the loading and unloading points are very small 
when the queues are not in existence. 

Figures 4 and 5 are of a slightly different type. Total time 
to unload ten thousand tons of cargo is again plotted on the vertical 
axis, as in Figures 2 and 3. The number of loading and unloading 
points is held constant at 4 each, their capacity is held constant at 
15 tons per hour, and the ship-to-shore distance is held constant at 
5 miles, as was done in Figures 2 and 3. 

In Figure 4, however, the number of transfer vehicles is held 
constant at 12, and vehicle payload is plotted against time. Four 
curves are shown, corresponding to vehicle speeds of 5, 10, 20, and 
40 knots. In the figure, S is used to denote speed. In Figure 5, 
the number of transfer vehicles is again held constant at 12 and 
vehicle speed is plotted against time. Four curves are shown, corre- 
sponding to vehicle payloads of 5, 10, 15, and 20 tons. In the figure, 
P is used to denote payload. The curves in Figures 4 and 5 display 
the same general shape as the curves in Figures 2 and 3. Again, the 
curves flatten out rapidly as we proceed from lower to higher values 
of the independent variables, reflecting the rapidly diminishing mar- 
ginal returns which are characteristic of this system. 
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FIGURE 4 



Time to Unload Ten Thousand Tons of Cargo versus Transfer 
Vehicle Payload (Distance 5 Miles, 12 Vehicles) 
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FIGURE 5 



Time to Unload Ten Thousand Tons of Cargo versus Transfer 
Vehicle Speed (Distance 5 Miles, 12 Vehicles) 
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All of the foregoing leads immediately to the realization that, 
in order to create an "efficient" system, or to improve on an existing 
system, we need to consider several variables simultaneously. The 
complexity of the relationships between the variables under discussion 
and the large number of possible combinations of parameter values that 
might have to be investigated lead to the conclusion that a simple 
function relating the variables would be of use even if it were only 
an approximating formula. Computational ease might compensate for 
lack of exactness, and once the approximating formula has led us to an 
"efficient" neighborhood, more detailed investigations can be carried 
out in this neighborhood by means of the computer simulation. 

Development of an Approximating Equation . 

In developing an approximating formula, it will be useful to 
start from a highly simplified model, from which a more sophisticated 
model may be developed. Let us consider a system in which the times 
required to load or unload a vehicle and the time for a vehicle to 
travel a given distance are known constants, rather than random var- 
iables. We will refer to this as a system with determinate times, or 
simply as a determinate system. Let us also assume that it is a bal- 
anced system, as discussed earlier, and that the fleet of transfer 
vehicles is composed entirely of vehicles of one type. 

For a balanced, determinate system the computation of total 
time to move a given tonnage of cargo through the system is fairly 
simple. We can readily design a balanced, determinate system in 
which there will be no queueing time. The total time to complete 
the movement of a given tonnage of cargo through such a system can 
then be determined as follows. Let T be the total tonnage of cargo 
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to be moved, let P be the payload of an individual transfer vehicle, 

D the ship--to-shore distance, S the speed of a transfer vehicle, N 
the number of transfer vehicles, and L the capacity of an individual 
loading (or unloading) point. Then the number of vehicle round trips 
required is the tonnage divided by the vehicle payload, or T/P. The 
number of round trips per vehicle is then T/NP. The one-way travel 
time is the ship-to-shore distance, D, divided by vehicle speed, S. 

The time required for a vehicle to make one round trip will be twice 
the one-way travel time plus the loading time plus the unloading time. 
Loading time and unloading time are assumed to be equal and are given 
by vehicle payload, P, divided by L, the capacity (in tons per hour) 
of an individual loading or unloading point. Thus the time required 
for a vehicle to make one round trip is 2 (D/S + P/L) . Multiplying 
round trip time by the number of round trips required per vehicle 
gives us approximately the total time required. 

Time = 2 (T/NP) (D/S + P/L) . (1) 

The formula is approximate in this case only because of the transient 
period at the start and at the end of system operation. 

Equation (1) assumed that we had designed a system in which 
the transfer vehicles never had to wait to load or to unload. A queue 
can develop in a determinate system with homogeneous vehicles only if 
the time required for a vehicle to travel to the unloading area, un- 
load, and travel back to the loading area is less than the time required 
to load the remaining N-1 vehicles. The time required for a vehicle to 
travel to the unloading area, unload, and travel back to the loading 
area is given by 2D/S + P/L. The time required to load the remaining 
N-1 vehicles is (N-1)(P/L) if we have one loading point. Only one queue 
will form in a balanced, determinate system, and since we started all 
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the transfer vehicles at the loading area, the one queue will be at 
the loading area <» Q, the amount of queueing time per vehicle round- 
trip, is given by 

Q = (N-1)(P/L) » (2D/S -h P/L). (2) 

If Equation (2) gives a negative value for Q, queueing time per 
vehicle round trip will be zero. 

Now consider a balanced determinate system with multiple load- 
ing and unloading points. Let K represent the number of loading points. 
Since we are dealing with a balanced system, K also represents the num- 
ber of unloading points. We will assume that all loading points oper- 
ate until the ship as unloaded, and that transfer vehicles depart and 
arrive in groups of K vehicles. A queue can develop only if the first 
group of K vehicles can travel to the unloading area, unload, and re- 
turn to the loading area before the remaining N-K vehicles have been 
loaded. The time required to load the remaining N-K vehicles will be 
given by ( (N-K) /K) (P/L) , and an approximate equation for queueing time 
per vehicle round trip in this system will be 

Q - ■ (P/L) - (2D/S + P/L). (3) 

By multiplying and grouping terms, we can rewrite this equation in the 
more convenient form 

Q = - 2 (D/S + P/L). 

We are now ready to consider the more realistic case of a sto- 
chastic model, in which the times for loading, unloading, and travel 
are random variables rather than constants. The amount of queueing 
time per vehicle round trip also becomes a random variable. We wish 
to develop a simple approximating equation for determining the average 
queueing time per vehicle round trip, which we shall denote by X. It 
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seems reasonable to suppose that the average queueing time per vehicle 
round trip for the stochastic model is related in some way to the 
queueing time per vehicle round trip for the determinate case. Values 
for Q, the amount of queueing time per vehicle round trip for the deter- 
minate model, were computed for a number of combinations of parameter 
values. Equation (3) was used for these computations. Values for X 
were then determined for the same combinations of parameters by means 
of the computer simulation. A least squares curve was then fitted to 
the data points. The relationship thus derived between X and Q is 
given by 

X = 0.19Q^ + 0.69Q + 0.37. (4) 

The root mean square error was 0.0456. This approximation gives 
reasonably accurate estimates of X for values of Q in the range from 
- 0.5 hours to + 1.5 hours. A negative value for Q indicates that in 
a determinate system there would be no queue and the loading and un- 
loading points would be unoccupied some of the time. Thirty data 
points were used in fitting the least squares curve. Table II shows 
the parameter values used, the average queueing time per vehicle round 
trip as determined by the computer simulation, and the average queue- 
ing time per round trip as determined by Equation (4). 

If the value of Q computed by Equation (3) is less than - 0.5 
hours, this suggests that the loading and unloading facilities are 
standing idle much of the time. If the computed value of Q is greater 
than 1.5 hours, the transfer vehicles are spending an inordinate amount 
of time waiting to load or to unload. In either case, the overall 
system operation would probably be considered inefficient. Thus, the 
range in which we will be primarily interested when applying the equa- 
tion will be the range covered by Equation (4), namely - 0.5 to + 1.5. 
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TABLE II 



Average Queueing Time per Vehicle Round Trip: 
Comparison of Simulation Results and ^ 
Results Computed by Approximating Equation 



Mean 


Mean 


Number 


Time 


Time 


Loading 


Travel 


of 


from 


from 


Time 


Time 


Vehicles 


Simulation 


Equation 


.06 Hr . 


.33 Hr . 


8 


.08 Hr . 


.29 Hr . 


.06 


.33 


12 


.35 


.52 


.06 


.33 


16 


.67 


.78 


.12 


.33 


8 


.06 


.21 


.12 


.33 


12 


.26 


.38 


.12 


.33 


16 


.56 


.70 


.25 


.33 


8 


.03 


- 0 - 


.25 


.33 


12 


.15 


.26 


.25 


.33 


16 


.36 


.49 


.50 


.33 


16 


.15 


.16 


.06 


.67 


8 


.24 


.29 


.06 


.67 


12 


.81 


.80 


.06 


.67 


16 


1.45 


1.48 


.12 


.67 


8 


.18 


.21 


.12 


.67 


12 


.77 


.70 


.12 


.67 


16 


1.32 


1.35 


.25 


.67 


8 


.12 


- 0 - 


.25 


.67 


12 


.53 


.49 


.25 


.67 


16 


1.18 


1,08 


.50 


.67 


12 


.27 


.16 


.50 


.67 


16 


.74 


,62 


.06 


1.0 


8 


.36 


.29 


,06 


1.0 


12 


1.17 


1.12 


.12 


1.0 


8 


.37 


.21 


.12 


1.0 


12 


1.14 


1.00 


.12 


1.0 


16 


2.04 


2.17 


.25 


1.0 


8 


.30 


- 0 - 


.25 


1.0 


12 


.95 


.76 


.25 


1.0 


16 


1.96 


1.83 


.50 


1.0 


12 


.61 
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The average time for a vehicle to complete a round trip in the 
stochastic model will be the sum of average travel times, average load- 
ing and unloading times, and X, the average queueing time per vehicle 
round trip. Average time per round trip, then, will be given by 
2 (D/S + P/L) -f (X). 

Total time required to unload a given tonnage of cargo will be approxi- 
mately the time per vehicle round trip multiplied by the required num- 
ber of round trips per vehicle, or 

Total Time = 2 (D/S + P/L) (T/NP) + (X) (T/NP) , (5) 

where X is given by Equation (4) . 

Application of the Approximating Equation . 

An application of the approximating equation would be to estimate 
the performance of a given ship off-loading system. The need for such an 
estimate might easily arise in the course of operational planning. In 
the same context of operational planning, we may perhaps find ourselves 
able to vary only one or two of the parameters which are treated as vari- 
ables in this model. The approximating equation will provide us with a 
ready means of estimating how changes in these parameters will affect 
the time required for the general unloading. In a somewhat longer time 
frame, we may be faced with a choice between several proposed systems. 

We could use either the computer simulation itself or the approximating 
equation to obtain estimates of the performance of each of the alterna- 
tive systems. We would then be able to choose the system that provides 
the best performance in relation to its cost. 

We have outlined the use of the computer simulation or the approx- 
imating equation in estimating the performance of a given system or in 
choosing between a given group of alternative systems. Another situation 
with which we may be faced is the task of determining the best way in 
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which to modify an existing system. By use of the approximating form- 
ula or the computer simulation we can estimate the reduction in time to 
be expected from incremental changes in each of the variables. If we 
also knew the additional cost associated with the proposed change in 
each variable, we could then select which variable to modify in order 
to achieve the largest reduction in time for a given cost. 

Suppose we are given a system in which the parameter values are 
as follows: the number of transfer vehicles, N, = 12, the speed of the 
transfer vehicles, S, = 10, the payload of the transfer vehicles, P, 

= 5, the number of loading points, K, = 4, the capacity of the loading 
points, L, = 15, and the ship-to-shore distance, D, = 5. From the data 
generated by the computer simulation we obtain an estimate of 287 hours 
required to move ten thousand tons of cargo through this system. For 
simplicity, suppose we wish to consider two possible changes in the sys- 
tem: first, a change in the payload of the twelve vehicles from five to 
ten tons, and second, a change in the number of five-ton payload vehicles 
from twelve to sixteen* Let us further suppose that the additional costs 
associated with these two changes are equal. We wish to determine which 
of the changes to implement in order to achieve the largest reduction in 
time. Changing the number of vehicles in the system from twelve to six- 
teen reduces the estimated time for completion of the unloading to 229 
hours. Changing the payload of the twelve vehicles from five tons to 
ten tons reduces the estimated time to 215 hours. In this hypothetical 
situation we would choose to increase vehicle payload rather than adding 
more vehicles to the system. 

A far more complex problem arises if we wish to design a complete 
system. There does not appear to be any global optimum in the system we 
are modeling. Decreasing ship- to-shore distance or increasing any one 
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of the other variables always seems to result in at least some reduc- 
tion in the time required to land a given tonnage of cargo. We can 
strive to design a system to accomplish the unloading in a minimum 
amount of time for a given cost or to accomplish unloading at a mini- 
mum cost in a given time. 

For illustrative purposes, suppose we wish to minimize the 
cost to land a given tonnage of cargo in a fixed amount of time. It 
appears that the most practical procedure would be simply to enumerate 
the combinations of parameter values that produce the desired result 
in terms of time to complete the landing of a given tonnage of cargo. 
This would appear to be a sizable task. However, the number of trans- 
fer vehicles in the system and the number of loading and unloading 
points clearly must be integer values. The speed and payload of the 
transfer vehicles and the capacity of the loading points might be 
thought of as continuous variables, but in practice we need only con- 
cern ourselves with a reasonable number of values for these variables. 

We might, for instance, consider five-ton increments in transfer vehicle 
payload. Even so, we will have to examine many combinations of para- 
meter values. Fortunately, a computer can solve the approximating 
equation very rapidly and can readily be programmed to eliminate all 
combinations that do not meet our time criterion. Many combinations 
can be eliminated almost automatically. For instance, if we find a 
set of parameter values including a fleet of twelve vehicles that will 
meet our unloading time criterion, we would not want to consider sys- 
tems using larger numbers of transfer vehicles in conjunction with the 
same set of values for the other parameters. Rough cost data should 
further serve to reduce the number of alternatives to manageable pro- 
portions. For this reduced number of alternatives we can then make 
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more detailed cost estimates and select the system that will accomplish 
the task at a minimum cost. 

We have developed an approximating equation for estimating the 
time required to land a given tonnage of cargo, and we have outlined 
how this equation or the computer simulation may be used in operational 
planning, in choosing between proposed alternative systems, and in deter- 
mining how best to modify an existing system. We have also described a 
procedure which offers a possible approach to the complex problem of 
total system design. A number of questions still remain, however. The 
data for the analysis in this chapter was based on the assumption of 
the exponential distribution for both service and travel time. The 
question of how this assumption affects system performance will be 
examined in Chapter V. In Chapter VI, we will examine some of the 
implications of having a fleet of transfer vehicles composed of more 
than one type of vehicle. 
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CHAPTER V 



TESTING THE CHARACTERISTICS OF THE SIMULATION MODEL 
A considerable amount of theoretical and practical work has 

g 

been done on networks of queues and on queues in tandem. This work 
generally assumes the existence of an independent stochastic process 
producing inputs or arrivals to the queueing network. Frequently, it 
is assumed that these arrivals occur in accordance with a Poisson dis- 
tribution. Because of the circular nature of the physical system which 
we are attempting to model, however, there is no independent input into 
the system. Rather, once the system is in operation the arrival rate 
at any point in the system is dependent on the output of the rest of 
the system. In order to gain some insight into the operation of the 
system and to determine whether changes in the underlying distributions 
of travel time and service time affect the operation, we will test two 
hypotheses about the system. First, we will test the hypothesis that 
the vehicles cycling through the system arrive at the loading area in 
accordance with a Poisson distribution. Secondly, we will test an 
hypothesis concerning the mean arrival rate. 

Testing Whether the System Is of the Poisson Type . 

In order to test whether arrivals at the loading area occur in 
accordance with a Poisson distribution, a means for measuring the 
arrival times at the loading area was inserted into the computer pro- 
gram. We wish to measure the amount of time, T, required to observe 

n arrivals. We will denote the time from the start of measurement 
t h. 

until the i arrival as U. . If the arrivals have occurred in accord- 

1 

ance with a Poisson distribution, then the random variables . . . , 

are independent and uniformly distributed over the interval (0,T) . 
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Since the number of observations, n, varied from five hundred to two 

thousand, we may use the central limit theorem to assert that the sum, 

n 



S = 
n 



2 

i=l 






of n independent random variables, each uniformly distributed on the 

interval 0 to T, may be considered to be normally distributed with t 

2 

mean n(T/2) and variance n (T /12). At a 95 per cent level of signifi- 
cance, then, we could accept the hypothesis that the arrivals are in 
accordance with a Poisson distribution if the following inequalities 
are satisfied: 



n(T/2) - (1.96) 



\l 



n(T /12) 



S ^ n(T/2) + (1.96) 



N 



n(T /12) 



This hypothesis was tested with data from 208 simulation runs. Var- 
ious distributions were assumed for both travel and service times and 
various values were used for the number of vehicles, number of loading 
points, mean travel time, and mean loading time, as shown in Table III. 

In each of the 208 tests, the hypothesis that the arrivals 
occurred in accordance with a Poisson distribution was accepted at a 
95 per cent level of significance. From these tests we conclude that 
the transfer vehicles are arriving at the loading area in accordance 
with a Poisson distribution. Furthermore, there is no indication that 
the Poisson nature of the arrivals in this system is dependent upon 
the underlying distributions of travel time and service time. We have 
not examined a large number of distribution functions, but those we 
have examined include the rather extreme cases of the uniform distribu- 
tion on the one hand and a constant value on the other. Although no 
skewed distributions were tested, it appears that the assumption that 
the arrivals occur in accordance with a Poisson distribution is a 
robust assumption. 
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TABLE III 



Distributions and Parameter Values 
Employed in Hypothesis Testing 



Distribution of 


Exp. 


Exp . 


Exp . 


Exp. 


Uniform 


Service Time 












Distribution of 


Exp . 


Exp . 


Constant 


Uniform 


Exp . 


Travel Time 












Mean 


0.33 Hr. 




0.33 Hr. 


0.33 Hr. 


0.25 Hr. 


Loading 


0.67 Hr. 


0.50 Hr . 


0.67 Hr. 


0.67 Hr. 


0.50 Hr. 


Time 


1.00 Hr. 




1.00 Hr. 


1.00 Hr . 


0.75 Hr. 




1.33 Hr. 








1.00 Hr. 


Mean 


0.06 Hr. 




0.12 Hr. 




0.12 Hr. 


Travel 


0.12 Hr. 




0.25 Hr. 




0.25 Hr. 


Time 


0.25 Hr. 




0.50 Hr. 




0.50 Hr. 




0.50 Hr. 


1.00 Hr. 


1.00 Hr. 


0.50 Hr. 


1.00 Hr. 




1.00 Hr. 




2.00 Hr. 




2.00 Hr. 




2.00 Hr. 












4.00 Hr. 










Number of 
Loading Points 


4 


1 

2 


4 


4 


4 


Number of 
Vehicles 


8,12,16 


2,3,4 

5,6 


8,12,16 


8,12, 16 


8,12, 16 


Number of 
Simulation Runs 


84 


10 


45 


9 


60 
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Testing an Hypothesis on the Mean Arrival Rate. 



Having accepted the hypothesis that arrivals at the loading 
area are occurring in accordance with a Poisson distribution, it next 
seemed desirable to make some determination concerning the mean value 
of this process. Intuitively, the mean arrival rate at the loading 
area should be the number of vehicles in the system divided by the 
average round-trip time, where the average round-trip time is taken 
to be the sum of the mean loading time, the mean unloading time, twice 
the mean travel time, the mean waiting time in the queue at the load- 
ing area, and the mean waiting time at the unloading area. If a 
Poisson process is observed until a predetermined number of events 
(in this case arrivals) have occurred, then the amount of observation 
time, W^, required to observe the n events can be used to test hypoth- 
eses concerning the mean value of the process. We shall denote the 
mean value by M and the hypothesized mean value by In order to 

avoid the transient period at the start of system operation, the mea- 
surement was made for 60 arrivals beginning at arrival number 400. 

If the events being observed are of the Poisson type, then 2MW is 

n 

chi-square distributed with 2n degrees of freedom. We shall accept 
the hypothesis that the mean value of the process is the hypothesized 
value Mq if the following inequalities are satisfied: 



'X ot/3 (2n) ^ ^ "X I - (2n) . 

2Mo n 2M„ 

This test was made on data from the same 208 simulation runs as 
described previously (Table III), and the hypothesis was accepted 
at a 95 per cent significance level on all but one of the 208 runs. 
(The one case in which the hypothesis was rejected appears to have 
been a statistical quirk, since it did not occur at an extreme value 
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or under any unusual conditions.) We conclude that the vehicles 
cycling through the system are arriving at the loading area at a mean 
rate equal to the number of vehicles in the system divided by the aver- 
age round-trip time, where the average round-trip time is the sum of 
the mean loading time, the mean unloading time, twice the mean travel 
time, the mean waiting time in the loading queue, and the mean waiting 
time in the unloading queue. Furthermore, there is no evidence from 
the simulation data that this mean arrival rate is affected by the 
underlying distribution of travel time or service time. 

Significance . 

We have accepted the hypothesis that, in the simulation model, 
arrivals at the loading area occur in accordance with a Poisson distri- 
bution, and we have accepted the hypothesis that the mean arrival rate 
at the loading area is determined by the number of vehicles in the sys- 
tem divided by the average round-trip time for a vehicle. We have 
tested these hypotheses over the range of values shown in Table III. 
There is no indication that these results are dependent on the assumed 
distribution of the travel times or the service times. Although these 
results throw some light on the nature of the process being investi- 
gated, they do not, of themselves, lead to an immediate closed form 
solution. The mean arrival rate, which is the parameter of the Poisson 
distribution, is not independent of, but rather depends upon the wait- 
ing times in the system. The waiting times are themselves dependent 
upon the arrival rate. Thus, if we knew either the arrival rate or 
the waiting times we could readily solve the system with closed form 
queueing equations and determine the values of the other. Before the 
fact, however, we do not know what the values will be for either the 
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mean arrival rate or the waiting times. Use of the equations developed 
in Chapter IV, however, will enable us to estimate the average queue- 
ing time per vehicle round trip and the total time required to complete 
the landing of a given tonnage of cargo. 

Up to this point we have been assuming that the fleet of trans- 
fer vehicles is composed entirely of one type of vehicle. Some of the 
complexities that arise in the case of mixed vehicle fleets will be 
examined in the next chapter . 
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CHAPTER VI 



IMPLICATIONS OF MIXED TRANSFER VEHICLE FLEETS 

We have on several occasions mentioned the fact that, in actual 
operations, our fleet of transfer vehicles consists of a variety of 
vehicle types rather than just one type. This is certainly true today, 
and seems likely to be the case in the future, due to the tactical 
problems that must be solved in the assault phase of any amphibious 
operation. Since this is so, it appears desirable to investigate the 
mixed transfer vehicle case, and our simulation provides a ready means 
for doing so. 

Fleets Composed of Vehicles with Differing Payloads . 

One of the situations that is of considerable interest is the 
case of a fleet of transfer vehicles of the same speed, but with differ- 
ing payloads. This is the situation that exists between various types 
of present-day landing craft (LCD’s, LCM’s, LCVP’s). The speed differ- 
ential is not significant, but there is a significant difference in the 
payloads of the various landing craft mentioned. 

We will consider two fleets of vehicles all having a speed of 
ten knots. In one fleet, half of the vehicles will have a fifteen-ton 
payload and the remainder a five- ton payload. In the other fleet, half 
of the vehicles will have a fifteen-ton payload and the remainder a 
ten-ton payload. When a number of vehicles are in the queue waiting 
to load or to unload, we will have a choice as to what type of vehicle 
to load or to unload first. One of our primary objectives in this 
portion of the study will be to determine an optimal decision rule for 
making this choice. Three possible decision rules which will be inves- 
tigated are: 
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1) load and unload the smaller payload vehicles first, 

2) load and unload the larger payload vehicles first, and 

3) choose vehicles in some random fashion for loading and 
unloading . 

In order to test these decision rules, the basic simulation pro- 
gram was modified slightly to include the decision rule as part of the 
simulation. In the basic simulation program, whenever a vehicle is to 
be brought from a queue into a loading or unloading point, the program 
searches through an array in which the staus of all vehicles is stored 
until it finds a vehicle whose status number indicates that that vehicle 
is in the loading or the unloading queue, as appropriate. Each vehicle 
has an index number in the sinulation. Among the vehicles in the queue, 
the simulation program always selects the vehicle with the lowest index 
number first. If the decision rule being tested was to load the smaller 
payload vehicles first, the payload assigned to the vehicles with the 
lower index numbers was the smaller payload, and vice versa if the deci- 
sion rule being tested was to load the larger vehicles first. Whenever 
a vehicle arrives and goes directly into a loading or an unloading point, 
its index number ajid hence its payload is known. When the decision rule 
being tested was to select vehicles according to some random procedure, 
even numbered vehicles were assigned one payload and odd numbered vehicles 
were assigned the other payload. 

In comparing the results achieved under the different decision 

rules, non-parametric or distribution free statistical tests were used. 

9 

The sign test and the signed rank test were selected. The hypothesis 
tested was the hypothesis of no difference between the results. Data 
were used only from combinations of parameter values that produced aver- 
age waiting times greater than 0.10 hour. When the average waiting time 
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is less than 0.10 hour per round trip, the decision rule is not playing 
an important role, since there is usually no queue or only one vehicle 
in the queue, and hence no choice exists as to which type of vehicle to 
load or to unload first. The data are shown in Table IV. 

In the case of the transfer vehicle fleet where half the 
vehicles had a fifteen-ton payload and half had a five-ton payload, the 
results obtained using the decision rule to load or unload the smaller 
vehicle first were tested against the other two decision rules described 
above. Using the signed rank test, the hypothesis of no difference in 
the results was rejected in both cases at the five per cent significance 
level. For the case of the fleet composed half of vehicles with a fifteen- 
ton payload and half of vehicles with a ten-ton payload, the results 
showed generally smaller differences, but again the hypothesis of no 
difference was rejected in both cases at the five per cent significance 
level . 

Other things being equal, larger payload vehicles (employed as 
pure fleets) are more efficient than smaller payload vehicles in that 
they move a fixed tonnage of cargo in a shorter period of time, as we 
have seen in Chapter IV. In considering the mixed fleet of vehicles, 
however, we conclude that we can improve the operation of the system 
by using the decision rule to load the smaller vehicles first when there 
is a choice, rather than loading and unloading the larger vehicles first 
or choosing vehicles at random. This improvement is apparently due to 
better utilization of the loading and unloading points. It should be 
borne in mind, however, that no allowance was made for buffer time in 
the simulation. In comparing the decision rule as applied to the fleet 
with fifteen-ton and five-ton vehicles, using the sign test, we would 
reject the hypothesis that loading and unloading the smaller vehicles 
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TABLE IV 



Time Required to Unload Ten Thousand Tons of Cargo with a Mixed Fleet 
of Transfer Vehicles: Comparison of Decision Rules When Half the 
Vehicles Have the Larger Payload and Half Have the Smaller 



Payload 

of 

Smaller 

Vehicle 


Mean 

Travel 

Time 


Number 

of 

Vehicles 




Choose 
15 Ton 
Vehic les 
First 


Choose 

Vehicles 

at 

Random 


Choose 

Smaller 

Vehicles 

First 


5 Tons 


0.12 Hr. 


12 




207 Hrs. 


197 Hrs. 


190 Hrs. 


5 


0.12 


16 




198 


193 


183 


5 


0.25 


12 




217 


209 


200 


5 


0.25 


16 




197 


207 


186 


5 


0.50 


12 




237 


232 


228 


5 


0.50 


16 




210 


211 


204 


5 


1.0 


16 




232 


244 


239 
















10 Tons 


0.12 Hr. 


8 




220 Hrs. 


216 Hrs. 


219 Hrs. 


10 


0.12 


12 




197 


193 


188 


10 


0.12 


16 




190 


190 


180 


10 


0.25 


8 




230 


242 


235 


10 


0.25 


12 




207 


198 


187 


10 


0.25 


16 




183 


188 


185 


10 


0.50 


12 




213 


210 


207 


10 


0.50 


16 




197 


191 


194 


10 


1.0 


12 




259 


253 


261 


10 


1.0 


16 




213 


214 


214 


10 


2.0 


16 




300 


299 


288 
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first produced a five per cent reduction in time compared to the deci- 
sion rule to load the larger vehicles first, but we would accept the 
hypothesis of a three per cent reduction. In percentage terms, the 
improvement is small, but the cost of implementing the decision rule 
is negligible. 

Fleets Composed of Vehicles with Differing Speeds . 

Another case we wish to consider is that of a mixed fleet in 
which the transfer vehicles all have the same payload but have differ- 
ing speeds. The procedure will be generally the same as that described 
above. The parameter values used and the results obtained from the simu- 
lation runs are shown in Table V. Using the sign test, we would accept 
they hypothesis that there is no difference between the decision rules 
at the five per cent significance level but we would reject the hypoth- 
esis at the ten per cent significance level. If the two decision rules 
we have considered do make any difference in the overall performance of 
the system under these conditions, this difference is apparently on the 
order of one per cent or less. Thus, we conclude that if the fleet of 
transfer vehicles consists of vehicles with equal payloads but differ- 
ent speeds, it does not make much difference whether we load the faster 
vehicles first or the slower vehicles first. 

Interpretation and Significance of the Results . 

In employing mixed fleets of transfer vehicles, it is apparent 
that part of the doctrine for their employment (the decision rules con- 
sidered here) can play an important role in the overall system perfor- 
mance in terms of time required to complete unloading. We have examined 
in detail just two of the numerous situations that might arise in prac- 
tice. However, for any given situation the computer simulation can 
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TABLE V 



Time Required to Unload Ten Thousand Tons of Cargo with a Mixed Fleet 
of Transfer Vehicles: Comparison of Decision Rules When 
Half the Vehicles Have a Speed of Ten Knots 
and Half Have a Speed of Twenty Knots 



Mean 

Loading 

Time 


Mean 

Travel for 
S lower 
Vehicles 


Number 

of 

Vehicles 


Choose 
20 Knot 
Vehicles 
First 


Choose 
10 Knot 
Vehicles 
First 


0.33 Hr. 


0.12 Hr. 


12 


188 Hrs. 


191 Hrs. 


0.33 


0.12 


16 


184 


185 


0.33 


0.25 


12 


202 


206 


0.33 


0.25 


16 


185 


186 


0.33 


0.50 


16 


200 


210 


0.67 


0.12 


8 


220 


221 


0.67 


0.12 


12 


194 


184 


0.67 


0.12 


16 


170 


187 


0.67 


0.25 


8 


232 


229 


0.67 


0.25 


12 


195 


194 


0.67 


0.25 


16 


177 


184 


0.67 


0.50 


12 


208 


212 


0.67 


0.50 


16 


179 


183 


0.67 


1.0 


12 


237 


245 


0.67 


1.0 


16 


199 


212 


0.67 


2.0 


16 


263 


268 
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readily be modified to reflect that situation and to test the results 
obtained by using alternative decision rules. The results are not 
obvious or easily obtainable by other means. Use of the simulation 
to determine the optimal decision rule appears to be a necessary part 
of the analysis whenever we are faced with the problem of attempting 
to optimize or to improve upon an existing system which includes a 
mixed fleet of transfer vehicles. It would be hazardous to attempt 
to extrapolate the results obtained outside the range of values actu- 
ally tested, but the simulation provides a ready means for testing 
various decision rules applied to any actual combination of vehicles. 

In dealing with the larger problem of designing a ship off- 
loading system, additional complexities are introduced by considera- 
tions of mixed fleets of transfer vehicles. From a cost effectiveness 
point of view, an optimal system might well be a system using a mixed 
fleet of transfer vehicles. By optimal we mean a system that would 
achieve a minimum unloading time for a given cost or a minimum cost 
for a given unloading time. To illustrate why a mixed fleet may be 
optimal, we will consider three cases, 

1) a pure fleet composed of twelve vehicles with five- 
ton payloads, 

2) a pure fleet composed of twelve vehicles with fifteen- 
ton payloads, and 

3) a mixed fleet composed of six vehicles with five-ton 
payloads and six vehicles with fifteen-ton payloads. 

The mixed fleet will be operated according to the optimal decision 
rule of loading and unloading the smaller payload vehicles first when- 
ever there is a choice. The data resulting from these tests are dis- 
played in Table VI. Fifteen different combinations of parameter values 
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TABLE VI 



Time Required to Unload Ten Thousand Tons of Cargo: Comparison of 
Simulation Results Achieved by Two Pure Fleets and a Mixed Fleet 



Mean 

Travel 

Time 


Number 

of 

Vehicles 


Pure Fleet 
of 

5 Ton 
Vehicles 


Mixed Fleet, 
Half 5 Ton 
Vehicles 
Half 15 Ton 
Vehicles 


Pure Fleet 
of 

15 Ton 
Vehicles 


0.12 Hr. 


8 


243 Hrs. 


230 Hrs. 


233 Hrs. 


0.12 


12 


196 


190 


186 


0.12 


16 


185 


183 


176 


0.25 


8 


292 


260 


235 


0.25 


12 


218 


200 


189 


0.25 


16 


188 


186 


185 


0.50 


8 


417 


317 


261 


0.50 


12 


287 


228 


197 


0.50 


16 


228 


204 


182 


1.0 


8 


661 


460 


334 


1.0 


12 


444 


289 


237 


1.0 


16 


346 


239 


200 


2.0 


8 


1154 


700 


491 


2.0 


12 


770 


445 


341 


2.0 


16 


592 


360 


263 
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were considered. The pure fleet of fifteen-ton vehicles achieved an 
average reduction in unloading time of 29.6 per cent as compared to 
the pure fleet of five-ton vehicles. The mixed fleet, however, did 
almost as well, achieving an average reduction in time of 20.5 per cent 
as compared to the pure fleet of five-ton vehicles. A mixed fleet of 
twelve vehicles, consisting of six five-ton vehicles and six fifteen- 
ton vehicles, was clearly superior in these simulation runs to a pure 
fleet of eight fifteen-ton vehicles but inferior to a pure fleet of 
twelve fifteen-ton vehicles. The twelve vehicle mixed fleet would be 
equivalent to a pure fleet of nine, ten or eleven of the fifteen-ton 
vehicles. Depending on the cost functions involved, it might well be 
that the mixed fleet would be the most cost-effective. This possibil- 
ity introduces considerable additional complexity to the already complex 
problem of designing a general unloading system. It appears at this 
point that a reasonable approach would be to use the procedure out- 
lined in Chapter IV to arrive at a cost-effective system employing a 
pure fleet, then test the implications of switching to a mixed fleet 
of transfer vehicles. 
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CHAPTER VII 



SUMMARY, SUGGESTED EXTENSIONS, AND GENERAL APPLICABILITY 

The design of an efficient system for the general unloading 
phase of an amphibious operation involves complex inter-relationships 
among a large number of variables. We have examined six of these var- 
iables, namely transfer vehicle speed and payload, the number of trans- 
fer vehicles in the system, ship- to-shore distance, and the number and 
capacity of loading and unloading points. The queueing problems that 
arise in the operation of this system play a highly significant role 
in overall system performance and hence must be taken into account. 

Due to the circular nature of the system we are attempting to 
model, closed form solutions are not readily attainable. We have struc- 
tured a model and developed a computer simulation by means of which the 
general unloading phase of an amphibious operation may be studied. 
Utilizing data generated by the computer simulation, we have examined 
some of the relationships between the variables in the system and we 
have developed an approximating equation for the time required to com- 
plete the unloading of a given tonnage of cargo. We have also investi- 
gated the impact of various hypothetical distributions of service time 
and travel time. The Poisson distribution of the arrivals does not 
appear to be dependent upon the underlying distributions of service 
time or travel time. The mean arrival rate, which is the parameter of 
the Poisson distribution, depends upon the mean values of service time 
and travel time, but does not appear to be dependent upon the probability 
distribution of service time or travel time. We have also examined the 
system when mixed fleets of transfer vehicles are employed. For two 
cases likely to arise in practice, we have compared the results obtained 
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using different decision rules as to which vehicles to load or unload 
first . 

There are several different levels at which we may be interested 
in applying the results of this paper. For convenience, we will discuss 
them in terms of the short-run time period, the mid-range time period, 
and the long-run time period. In the short run we may be faced with 
an existing sytem, in which case our only recourse is to employ the 
most effective decision rule in operating the system. We have illus- 
trated how a system employing a mixed fleet of transfer vehicles can be 
analyzed in order to achieve the best possible results with fixed re- 
sources. In the context of operational planning, we have outlined how 
the approximating formula or the computer simulation itself may be used 
to obtain estimates of the performance of a given system. 

In applications pertaining to the mid-range time period from 
five to ten years hence, we may be faced with the task of determining 
the best means of improving upon an existing system or we may be faced 
with the task of choosing among a number of alternative systems. In 
Chapter IV we outlined procedures for dealing with these situations. 

In the long run we are free to modify all of the variables under 
consideration. This presents us with by far the most complicated prob- 
lem. We have defined an efficient system to be a system in which the 
ratio of marginal return (in time reduction) to marginal cost is the 
same for each of the variables. We have mentioned some of the practical 
and theoretical problems involved in trying to find such efficient 
points. Nevertheless, the step-by-step procedure we have outlined in 
Chapter IV appears to be the best practical approach to achieving an 
efficient system design. 
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Prior to implementing the results developed in this paper, it 
would be necessary to test the entire model by comparing the results 
predicted by the model with results achieved in actual operations. It 
would also be necessary to determine the actual parameter values to be 
used. This should be accomplished by field experimentation under real- 
istic operating conditions. It appears that accurate estimates of the 
parameter values for loading and unloading time and for travel time are 
more important than the probability distributions associated with those 
t ime s . 

Ship characteristics determine the number and capacity of the 
loading points in the system which we have modeled. Lead times for 
ship design and construction are such that the characteristics of the 
ships that will be available five to ten years in the future are fairly 
well known today. Thus we might reasonably be able to fix the number 
and capacity of loading points that will be available in a "typical" 
task force in the time period five to ten years hence. With appropriate 
cost data, and after validation of the model, it should be possible to 
develop transfer-vehicle fleets that meet our definition of an efficient 
system for those variables not considered fixed. Unless some unexpected 
dominant transfer vehicle fleet emerges, it would probably be necessary 
to make separate analyses for various ship- to-shore distances. The 
implications of various ship- to-shore distances on system cost and per- 
formance could then be portrayed. Alternatively, it might be possible 
to assign a discrete probability mass function to the variable ship-to- 
shore distance and' then to develop a transfer vehicle fleet that would 
be efficient in an expected value sense. That is, it would be efficient 
over a large number of operations if the various ship- to-shore distances 
occurred in accordance with the assigned probability mass function. 
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Another interesting extension would be to compare the per- 
formance of a system employing amphibian vehicles to move cargo from 
the ships to an inland dump and a system employing landing craft and 
trucks to accomplish the same task. The model and simulation we have 
developed can serve as the model of amphibian vehicles shuttling from 
a ship to an inland dump or as a model of landing craft shuttling from 
a ship to the waterline at the beach. A system employing landing craft 
and trucks to move cargo from a ship to an inland dump involves the 
interaction, at the beach, of two circular queueing systems such as the 
one we have modeled. A comparison of the relative effectiveness of a 
system employing amphibian vehicles and the system employing landing 
craft and trucks could then be made. This would be of considerable 
interest since a major tactical and logistical consideration in amphib- 
ious operations is to move men, equipment, and supplies rapidly inland 
and to avoid congestion on the beach itself. 

This model and the accompanying computer simulation were devel- 
oped in order to investigate the general unloading phase of an amphib- 
ious operation, but the model would be generally applicable to any 
closed system in which a fixed population cycles continually through 
the system. Many mater ials-handling systems and systems incorporating 
shuttling operations would appear to be of the same general nature and 
hence amenable to analysis by the same general methods and by the 
same sort of model and simulation. 
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APPENDIX 



The simulation program was written in FORTRAN IV for 
the IBM 360/67 computer. This appendix includes a list of the 
variable names used in the program, generalized flow charts, 
and a computer print-out of the program. 
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LIST OF VARIABLE NAMES 



AL 


Mean Loading Time 


ALN 


Transfer Vehicle Payload 


AM 


Mean Travel Time 


AQL 


Average Waiting Time per Vehicle Round Trip at Loading Area 


AQU 


Average Waiting Time per Vehicle Round Trip at Unloading 
Area 


D 


Total Number of Round Trips Required 


KK 


Number of Loading (Unloading) Points 


KQUE 


Length of the Queue at the Unloading Area 


KS(K) 


Status of Vehicle K 


KSTAT (L) 


Status of Unloading Point L 


KVEH 


Number of Transfer Vehicles 


LQUE 


Length of the Queue at the Loading Area 


LSTAT (L) 


Status of Loading Point L 


QLL(J) 


Length of Time Loading Queue Was of Length J 


QLU(J) 


Length of Time Unloading Queue Was of Length J 


TIM 


Time 


TNARL(K) 


Time of Next Arrival at Loading Area for Vehicle K 


TNARU(K) 


Time of Next Arrival at Unloading Area for Vehicle K 


TNDPL(K) 


Time Next Departure from Loading Area for Vehicle K 


TNDPU(K) 


Time Next Departure from Unloading Area for Vehicle K 


TLE 


Length of Time at Least One Loading Point Was Empty 


TUE 


Length of Time at Least One Unloading Point Was Empty 


X 


Clock for Unloading Queue 


XX 


Clock for Unloading Point 


Y 


Clock for Loading Queue 
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YY Clock for Loading Point 

Z Tons of Cargo Moved 



61 




62 






63 






64 
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FORTRAN PROGRAM LISTING 



C 



60 



62 



65 

63 

2 

71 



DIMENSION 0LLI50) .OLU(50) fTNOPLI 50),TNOPU(50) 
1TNARU(50) ,TNARL(50),KSTAT(10I ,LST4T(10) ,KS(50 
IX = 123456709 
CALL RANDU ( IX, IY,RN) 

IX = lY 

INITIALIZE PARAMETERS AND ARRAYS FOR NEXT RUN 
KK — A 

00 103 IK = 1,4 
B = IK 

ALN = 5. ♦ B 
AL = ALN / 15. 

DO 106 IJ = 1,7 

AM = HI. / 32. ) ♦ 2. ♦♦ IJ) 

DO 59 I = 8,16,4 
JJ = I ♦ 1 
DO 60 K = 1, JJ 
OLL(K) = 0. 

QLU(K) = 0. 

TNOPL(K) = 9999. 

TNOPU(K) = 9999. 

TNARL(K) = 9999. 

TNARIMK) = 9999. 

KS(K) = 1 
CONTINUE 
U = 0. 

SN = 0. 

UI = 0. 

TIM = 0. 

X = 0. 

Y = 0. 

XX = 0. 

YY = 0. 

TUE = 0. 

TLE = 0. 

KOUE = 0 
Z = 0. 

DO 62 K = 1, KK 

KSTAT(K) = 0 

LSTAT(K) = 0 

CONTINUE 

LOUE = I 

DO 63 K = 1, KK 

CALL RANDU (IX,IY,RN) 

IX = lY 

R = (-AL ♦ ALOG(RN)) 

TNOPL(K) = TIM + R 

KS(K) = 2 

LSTAT(K) = 1 

LQUE = LOUE - 1 

IFILOUE .EO. 0) GO TO 2 

CONTINUE 

IF(Z .GE. 10000) GO TO 90 
L = 1 

DO 71 LK = 2, I 

IF(TNDPL(L) .LE. TNOPL(LK)) GO TO 71 
L = LK 
CONTINUE 
M = 1 

DO 72 MK = 2, I 

IF(TNARU(M) .LE. TNARU(MK)) GO TO 72 



I 
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c 



c 



72 

73 

74 

3 

4 
10 



14 



12 

11 



20 



22 

23 



M * MK 
CONTINUE 
LJ » 1 

DO 73 LK = 2,1 

IFCTNDPU(LJ) .LE. TNOPUILK)) GO TO 73 
LJ * LK 
CONTINUE 
MJ = 1 

DO 74 MK * 2, I 

IFCTNARLIMJ) ,LE. TNARLIMKI) GO TO 74 

MJ = HK 

CONTINUE 

IFITNDPUL) .GT, TNARUfMII GO TO 3 
IFItNDPUL) .GT. TNOPUILjn GO TO 4 
IFfTNDPL(L) .GT. TNARLIMJH GO TO 30 
GO TO 40 

IF(TNARU(H) .GT. TNDPU(LJI) GO TO 4 
IF(TNARU(M) .GT. TNARL(MJ)) GO TO 30 
GO TO 10 

IFITNDPUILJ) .GT. TNARLCMJ)) GO TO 30 
GO TO 20 

ARRIVAL AT UNLOADING AREA IS NEXT EVENT 
TIM = TNARU(M) 

IFIKQUE .NE. 0) GO TO 11 
DO 12 K = 1,KK 
IF(KSTATIK) .NE. 0) GO TO 12 
CALL RANDU (IX,IY,RN) 

IX = lY 

R = <-AL ♦ ALOGIRN)) 

TNDPUIMI = TNARU(M) + R 
KSTATIK) = 1 
KS(M) » 5 

TUE = TUE + (TNARU(M) - XX) 

XX « TNARUIM) 

J « KQUE * 1 

OLU(J) = OLUU) ♦ (TNARUCM) - X) 

X » TNARU(M) 

TNARU(H) * 9999. 

GO TO 2 
CONTINUE 

J = KOUE ♦ 1 

QLUIJ) » OLUU) ♦ ITNARU(M) - X) 

X * TNARUIM) 

KQUE s kQUE * 1 
TNARU(M) * 9999. 

KS(M) e 4 
GO TO 2 

DEPARTURE FROM UNLOADING AREA IS NEXT EVENT 
TIM = TNDPUCLJ) 

CALL RANDU (IX,IY,RN) 

IX = lY 

R = C-AM* ALOGIRN)) 

TNARL(LJ) = TIM ♦ R 

ksTlj) = 6 

IFUQUE .EQ. 0) GO TO 21 
J « KQUE + 1 

OLUU) » QLUU) ITIM - X) 

KQUE » KQUE - I 

X * TIM 

DO 22 K = 1. I 

IFIKSIK) .EQ. 4) GO TO 23 

CONTINUE 

KS(K) * 5 

CALL RANDU (IXtlY.RN) 

IX « lY 
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R = (->SL * ALOG(RN)) 

?7 TNDPUCO = ri'l + R 

TNDPUC J) = 99Q9. 

Z = Z ♦ ALN 
GO TO 2 

21 XX = T\iJPU(LJ) 

00 2A K = 1,KK 

IF(KSTA,{K) ,E0. 0) GO TO 24 
KSTAT(K ) = 0 
GO TO 25 
■’4 CONTINUE 
25 ’ Z + ALN 

29 J KOUE + 1 

OLU(J) = QLU(J) + (TNOPUILJ) - X) 

X ••= TNOPU(LJ) 

TNDPU(LJ) = 99f^9. 

GO TO 2 

C ARRIVAL AT LOADING AREA IS Ni-XT EVENT 

30 TIM = TNARL(MJ) 

36 IFILQUE .NE. 0) GO TO 31 
DO 32 K = 1,KK 

IF(LSTAT(K) .NE, 0) GO TO 32 
CALL RANDU ( IX, IY,RN) 

IX = lY 

R = (-AL ♦ ALOG(RN)) 

30 TNDPL(MJ) = TIM ♦ R 
LSTAT(KJ = 1 
KS(MJ) = 2 

TLE = TLE + (TIM - YY) 

YY = TIM 
J = LQUE + 1 

OLL(J) = OLL(J) + (TIM - Y) 

Y = TIM 

TNARL(MJ) = 9999. 

GO TO 2 
32 CONTINUE 

31 J = LOUE + 1 

OLL(J) = OLL(J) ♦ (TIM - Y) 

Y = TIM 

LQUE = LOUE + 1 
TNARL(MJ) = 9999. 

KS(MJ) = 1 
GO TO 2 

C DEPARTURE FROM LOADING AREA IS NEXT EVENT 

40 TIM = TNDPL(L) 

CALL RANDU ( IX, lY ,RN) 

IX = lY 

R = (-AM* ALOG(RN) ) 

TNARU(L) = TIM + R 
TNDPL(L) = 999C. 

KS(L) = 3 

IF(LOUE .EO. 0) GO TO 41 
J = LOUE + 1 

OLL(J) = QLL(J) + (TIM - Y) 

LQUE = LOUE - 1 

Y = TIM 

DO 42 K = 1, I 

IF(KS(K) .EO. 1) GO TO 43 

42 CONTINUE 

43 KS(K) = 2 

CALL RANDU (IX,IY,RN) 

IX = lY 

P = (-AL * ALOG(RN)) 

47 TNDPL(K) = TIM + R 
GO TO 2 
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41 



44 

45 



90 

91 



YY = TIM 
00 44 K = 1 
IFILSTAKK) 
LSTAT (K) = 
GO TO .45 
CONTINUE 
J = LOUE 
QLUJ) = 

Y = TIM 
GO TO 2 
KVEH = I 
PRINT 91 



, KK 

.EQ. 

0 



0) GO TO 44 



♦ 1 

OLL( J) 



♦ (TIM - Y) 



TIM, 
(//lOX, 



— I 



61 



94 



110 



105 

92 

59 

106 
103 



,F10,5, 

0 . 

0 . 

/ ALN 
J = 1,1 
- 1 

♦ QLUJ) 

* QLU(J) 



/ D 
/ D 
AQL, 



KVEH, TL£,TUE 
' time = • »F10,5, 



• TUE = ' ,F10.5) 



AQL 

AQU 



* 

♦ 



FORMAT 
1* TLE 
AQL = 

AQU = 

0 = Z 
DO 61 
A = J 
AQL = 

AQU = 

CONTINUE 
AQL « AQL 
AQU = AQU 

PRINT 94, AQL, AQU 

FORMAT (//lOX, • AVG DELAY LOADING = • 

1 • AVG DELAY UNLOADING = • ,F15.5) 

PRINT 110,AM,AL 
F0RMAT(//10X, 'AM • ,F10.5, • AL • ,F10.5) 
JJ = I ♦ 1 

DO 105 J = 1, JJ 

OLUCJ) = QLU(JI / TIM 

OLL(J) = QLL(J) / TIM 

CONTINUE 

PRINT 92, ( (QLU J),QLU(J) ), 

FORMAT (//10X,2F10.5) 

CONTINUE 
CONTINUE 
CONTINUE 
RETURN 
END 



• NO.VEH = • ,15, 



,F15,5 



J = 1,JJ) 
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